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Abstract. We develop a new variational formulation of the inverse Stefan problem, where information on the 
heat flux on the fixed boundary is missing and must be found along with the temperature and free boundary. 
We employ optimal control framework, where boundary heat flux and free boundary are components of the 
control vector, and optimality criteria consists of the minimization of the sum of L2-norm declinations from 
the available measurement of the temperature flux on the fixed boundary and available information on the 
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r^ ' phase transition temperature on the free boundary. This approach allows one to tackle situations when the 

phase transition temperature is not known explicitly, and is available through measurement with possible 
error. It also allows for the development of iterative numerical methods of least computational cost due to 
the fact that for every given control vector, the parabolic PDE is solved in a fixed region instead of full 
free boundary problem. We prove well-posedness in Sobolev spaces framework and convergence of discrete 
optimal control problems to the original problem both with respect to cost functional and control. 
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1 Description of Main Results 

1.1 Introduction and Motivation 

Consider the general one-phase Stefan problem ([lU [25]): find the temperature function 
u{x,t) and the free boundary x = s{t) from the following conditions 

{a{x,t)ux)x + b{x,t)ux + c{x,t)u — Ut = f{x,t), for (x,t)Gi7 (1.1) 

u{x, 0) = (j){x), 0<x< s(0) = So (1.2) 

a{0,t)u,{0,t) = g{t), 0<t<T (1.3) 

ais{t),t)uMt),t) + 7(s(t), t)s'(t) = x{s{t),t), 0<t<T (1.4) 

u{s{t),t) = fi{t), 0<t<T (1.5) 

where a, b, c, f, 0, g, 7, X; A* are known functions and 

a{x,t) > ao> 0, So > (1.6) 

Q = {{x,t) : 0<x < s{t), 0<t<T} 

In the physical context, / characterizes the density of the sources, is the initial temperature, 
g is the heat fiux on the fixed boundary and n is the phase transition temperature. 

Assume now that some of the data is not available, or involves some measurement error. 
For example, assume that the heat fiux g{t) on the fixed boundary x = is not known and 
must be found along with the temperature u{x,t) and the free boundary s{t). In order to 
do that, some additional information is needed. Assume that this additional information is 
given in the form of the temperature measurement along the boundary x = 0: 

u{0,t) = u{t), forO<t<r (1.7) 

Inverse Stefan Problem (ISP): Find the functions u{x,t) and s(t) and the boundary heat 
flux g it) satisfying conditions / (i.i|) - pT?[ ). 

ISP is not well posed in the sense of Hadamard. If there is no coordination between 
the input data, the exact solution may not exist. Even if it exists, it might be not unique, 
and most importantly, there is no continuous dependence of the solution on the data. In- 
verse Stefan problem was first mentioned in [5], in the form of finding a heat fiux on the 
fixed boundary which provides a desired free boundary. This problem is similar to non- 
characteristic Cauchy problem for the heat equation. The variational approach for solving 
this ill-posed inverse Stefan problem was performed in [6l [7]. First result on the optimal 
control of the Stefan problem appeared in [35]. It consists of finding optimal value of the 
external temperature along the fixed boundary, in order to ensure that the solutions of the 
Stefan problem are close to the measurements taken at the final moment. In [35] existence 
result was proved. In [37] the Frechet differentiability and the convergence of the differ- 
ence schemes was proved for the same problem and Tikhonov regularization was suggested. 
Later development of the inverse Stefan problem was along these two lines: Inverse Stefan 
problems with given phase boundaries were considered in [H [21 [Bl [IHl [HI [121 [13 [211 [15] ; 
optimal control of Stefan problems, or equivalently inverse problems with unknown phase 
boundaries were investigated in^[T3l[IHl[l9l[20l[ni[22l[Ml[28l[26l[29l[30l[3ll[l5]. We 
refer to monography [15] for a complete list of references of both types of inverse Stefan 



problems, both for linear and quasilinear parabolic equations. The main methods used to 
solve inverse Stefan problem are based on variational formulation, method of quasi-solutions 
or Tikhonov regularization which takes into account ill-posedness in terms of the dependence 
of the solution on the inaccuracy involved in the measurement (11.71) , Frechet differentiability 
and iterative conjugate gradient methods for numerical solution. Despite its effectiveness, 
this approach has some deficiencies in many practical applications: 

• Solution of the inverse Stefan problem is not continuously dependent on the phase 
transition temperature /i(t): small perturbation of the phase transition temperature 
may imply significant change of the solution to the inverse Stefan problem. Accordingly, 
any regularization which equally takes into account instability with respect to both z/(t) 
from measurement fll.7p . and the phase transition temperature fi{t) from (11.51) will be 
preferred. It should be also mentioned that in many applications the phase transition 
temperature is not known explicitly. In many processes the melting temperature of 
pure material at a given external action depends on the process evolution. For example, 
gallium (Ga, atomic number 31) may remain in the liquid phase at temperatures well 
below its mean melting temperature ([25j). 

• Numerical implementation of the iterative gradient type methods within the existing 
approach requires to solve full free boundary problem at every step of the iteration, and 
accordingly requires quite high computational cost. Iterative gradient method which 
requires at every step solution of the boundary value problem in a fixed region would 
definitely be much more effective in terms of the computational cost. 

The main goal of this project is to develop a new variational approach based on the optimal 
control theory which is capable of addressing both of the mentioned issues and allows the 
inverse Stefan problem to be solved numerically with least computational cost by using 
conjugate gradient methods in Hilbert spaces. In this paper we prove the existence of the 
optimal control and convergence of the family of time-discretized optimal control problems 
to the continuous problem both with respect to cost functional and control. We employ 
Sobolev spaces framework which allows to reduce the reguarity and structural requirements 
on the data. We address the problems of convergence of the fully discretized family of 
optimal control problems, Frechet differentiability and iterative conjugate gradient methods 
in Hilbert spaces in an upcoming paper. 

Throughout the paper we use usual notation for Sobolev spaces according to references 

In the next section we formulate a new variational formulation of the inverse problem 
which takes into account the described deficiencies. 

1.2 Optimal Control Problem 

Consider a minimization of the cost functional 

J{v) = Po\HO,t) - z/(t)||i^[o,T] +/3ilk(^W,i) -/^WIlL[o,r] (1-8) 

on the control set 
Vr = {v = is,g) G W^[0,T] X W^[0,T] : S < s(t) < /,s(0) = So,max( \\s\\wf \\g\\wi < R] 

where 6, 1, R, f3o, f3i are given positive numbers, and u = u{x, t; v) be a solution of the Neu- 
mann problem fll.ip - fll.4l) . 



Definition 1.1 The function u G W2'' {^) is called a weak solution of the problem ( I -?.-?)) - 
( fTTl J ifu{x,0) = 0(s) e 1^2^0,50] and 



= / / [aMx*^x - bu^^ - cu^ + Ut^ + /$] rfx dt 
Jo Jo 

[^{s{t),t)s\t)-x{s{t),t)]^{s{t),t)dt+ [ gmiO,t)dt (1.9) 

^0 

for arbitrary $ G IV2 ' (^) 

We also need a notion of weak solution from V2(fi) of the Neumann problem: 

Definition 1.2 The function u G V2(fi) is called a weak solution of / li.ip - PT^ «/ 

0= / / [aui.$a; - 6mx$ -CM$ -M<l>t + /<l>]dxdt- / 0(x)$(x, 0) (ix+ 
io io io 

(7(t)$(0,t)rft+ / [^{s{t),t)s\t) ~ u{s{t),t)s'{t) - x{s{t),t)]^{s{t),t)dt (1.10) 
^/o 

for arbitrary $ G VF2' (^) snc/i t/iat $|j=y = 0. 

If u is a weak solution either from V2(r2) (or W2' (^)), then traces m|^=o ^^"^ '"L=s(t) ^^^ 
elements of L2[0,T], when s G W2^[0,T] ([271 [23]) and cost functional J{v) is well defined. 
Furthermore, formulated optimal control problem will be called Problem /. 

1.3 Discrete Optimal Control Problem 

Let 

^r = {tj = j -r, j = 0, l,...,n} 

be a grid on [0,T] and r = — . Consider a discretized control set 

Vn = {[v]n = {[sU [9]n) e M2-+2 : 0<S<s,<l, max(||H„||^.; ||[^].||^,) < R'} 

where, 

[s]n = (so,si,...,s„) G M"+\ [g]n = {go,gi,...,gn) e M"+i 

n— 1 n n— 1 n— 1 n 

IIMnll^l = ^ r4 + ^ rsf, + 5^ r4 ,„ || [g]n\\l. =J2rgl + J2 ^^ffc- 

fe=0 fc=l fc=l fc=0 fc=l 

under the standard notation for the finite differences: 

_ Sk — Sfc„i _ Sfc+l — Sk 2 _ -^fc+l ^ ^-^fe + ^k-1 

St,k — - ) St,k — - , %^fc - ■;^ • 

Introduce two mappings Q„ and Vn between continuous and discrete control sets: 

Qn{v) = [v]n = {[s]n, [g]n) , for f G Vr 

where Sk = s{tk),gk = g{tk), A; = 0, 1, ..., n. 

Vn{\v]n) = t'" = {s\ g-) G iy|[0, T] X Wl[Q, T] for [v]^ G V^, 
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where 



Sk-i + {t-tk-i-l)sju^ + \{t-tk-ifsitui tk-i<t <tk,k = 2,n. ■ ' 



dk — dk-i 



5'"W = 9k-i H {t - tk-i), tk-i <t<tk,k = l,n. 

T 

Introduce Steklov averages 

1 Z"** 1 /■*'= 

dk{x) = — / d{x,t)dt, hk = — h(t)dt, 



^ Jtk-i ^ Jtk-i 



where d stands for any of the functions a, b, c, /, and h stands for any of the functions u, fi. 
Given v = {s,g) E Vr we define Steklov averages of traces 

xj = - r x{s{t),t)dt, i^ssT = - r i{s{t),t)s'{t)dt. (1.12) 



Given [v]n = {[s]n, [g]n) £ V^ we define Steklov averages Xs" ^^^ (7s" (s")')*^ through 01.12p 
with s replaced by s" from (II. lip . 

Next we define a discrete state vector through time-discretization of the integral identity 
(ILSD 



Definition 1.3 Given discrete control vector [v]n, the vector function 

[u{[v]n)]n = {u{x; 0), u{x; 1), ..., m(x; n)) 
is called a discrete state vector if 
(a)u{x;0) = (P{x)eWi[0,So]; 

(b) For arbitrary k = 1,2, . . . ,n, u{x; k) G W2^[0, Sk] satisfy the integral identity 

\ak{x) — Y — ^\^) ~ ^k — T^ — ri{x) - Ck{x)u{x; k)r]{x) + fk{x)r]{x) 

+ uj{x; k)r]{x)'^ dx + ((7s"(s")')' - x'")^(sfc) + ^^//(O) = 0, (1.13) 

for arbitrary rj G 1^2^ [0, Sk], where 

u(x; k) — u(x; k — 1) 
Ut{x; k) = . 

(c) For arbitrary k = 0,1, ..., n, u{x; k) G W2^[0, Sk] iteratively continued to [0, 1] as 

u{x; k) = u{Tsk - x;k), 2"~^Sfc < x < 2"sfc, n = l,nk, Uk < N = 1 + loga 7 (1-14) 
where [r] means integer part of the real number r. 



Consider a discrete optimal control problem of minimization of the cost functional 

n 2 " 2 

Xn(M„) = /3or5^(u(0;/c)-z/fe) + ftr J^ (u(sfc; A:) - /i^) (1.15) 



k=l k=l 



on a set V^ subject to the state vector defined in Definition 1.3. Furthermore, formulated 
discrete optimal control problem will be called Problem /„. 

Throughout we use piecewise constant and piecewise linear interpolations of the discrete 
state vector: given discrete state vector [n([f]„)]„ = (m(x; 0), m(x; 1), ..., u{x; n)), let 



u^(x, t) = u{x; k), if tk-i < t < tk, < x < I, k = 0,n, 



u'^{x,t) = u{x; A; — 1) + Uj{x; k){t — tk-i), if t^-i < t < tk, < x < I, k = l,n, 
u^{x, t) = u{x; n), ii t >T, < x < I. 
Obviously, we have 

1.4 Formulation of the Main Result 

Let 

D = {{x,t) : 0<x</, 0<t<r} 

Throughout the whole paper, with the exeption of Section 13.11 we assume the following 
conditions are satisfied by the data: 

a,b,ceLooiD), f eL2{D), 

(l)eW^[0,So], ^,x^W^^\D), /i,z/GL2[0,T], 

the coefficient a satisfies (11. 6p almost everywhere on D, the generalized derivatives %, || 
exists and 



T 



da 



dt 



dt < +00. (1.16) 



— G Loo(-D), / esssupQ<x<i 

Our main theorems read: 

Theorem 1.1 The Problem I has a solution, i.e. 

V, = {veVR: J{v) = J,= inf J{v)} ^ 

v&Vr 

Theorem 1.2 Sequence of discrete optimal control problems In approximates the optimal 
control problem I with respect to functional, i.e. 

lim X„, = X, (1.17) 

n— >+oo 

where 

In, = inf X„([t;]„), n = 1, 2, ... 

If [v]n^ E VJl is chosen such that 

^n, < X„([f]„J < In, + en, €„ | 0, 



then the sequence Vn = {sn,gn) = 'Pn{[v]nJ converges to some element v^ = {s^:,g^:) G K 
weakly in W2[0, T] x 14^2^ [0, T], and strongly in VF^j^fO, T] x L2[0, T]. In particular s„ converges 
to s* uniformly on [0,T]. Moreover, piecewise linear interpolation vT' of the discrete state 
vector [u[v]n^n converges to the solution u{x,t]v^) e W2' (fi*) of the Neumann problem 
(EjP-llip weakly m Wl'^iSl^). 



2 Preliminary Results 

In a Lenima [?!T] below we prove existence and uniqueness of the discrete state vector [M([f]n)]n 
(see Definition 11.31) for arbitrary discrete control vector [f ]„ G V^. In a Lemma 12.21 we 
remind a general approximation criteria for the optimal control problems from ([36]). In a 
Lemma 12.31 we prove some properties of the mappings Q„ and Vn between continuous and 
discrete control sets. 

Lemma 2.1 For sufficiently small time step r, there exists a unique discrete state vector 
[u{[v]n)]n for arbitrary discrete control vector [f]„ G V^. 

Proof. To prove uniqueness, it is enough to show that if 

u{x; k-l) = 0, (7."(s")')' = 0, x'n =0, g'' = 0, fk{x) = 

then u{x;k) which solves (I1.13P vanishes identically. Under these assumptions by choosing 

r]{x) = u{x; k) in (I1.13P we have 

'' fafc(a;) f ^^I^Ml)' _ ^^ ^Ma:;A:) ^ _ ^^^^^^2^^. ^^ ^ 1^2(^. ^)^ ^^ ^ q_ ^2.1) 

iQ \ \ dx J dx r / 

Using (II. 6p and Cauchy inequality with e > we derive that 

r" /du{x;k)\^ ^ 1 P 2/ ,Nj / 

Qq / ; 1 dx -\ — / u [x; k) dx < 

Jo \ dx J T Jq 

eM r" fdu{x;k)\2 ^m ,\ T* 2, ,,, , , 

where 

M = max [\\a\\L^^Dy, \\b\\L^iD); ||c||l^(d)J. 

By choosing e = a^/M in (12.20 we have 

2 Jq \ dx J \T TqJ Jq 

where 

^" = (2^ + ^'^ • 

From (12. 3 p it follows that u(x; /c) = if r < tq. 

To prove an existence we apply Galerkin method. Consider an approximate solution 

N 

un{x) = y^^djipi^x) 
1=1 
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where {'ipi} is a fundamental system in 1^2^ [0, Sk] and the coefficients {di} solve the following 
system 



Sk 



ak{x)——%l)[{x) - hk{x)-^^i{x) - Ck{x)uN{x)iJi{x) + -UN{x)ipi{x) + fk{x)tpi{x) 

1 



dx 



'' u{x; k - l)Ux) dx - [(7."(s")')' - X'n] Usk) - 9%{0) t = l,...,N (2.4) 



which is equivalent to 



y^ / ak{x)i)j{x)i)-{x) -bk{x)i)j{x)i!i{x) - Ck{x)i)j{x)'ipi{x) + -i)j{x)i!i{x) dx dj 



i=i 



Sk 



fk{x)Mx) + -u{x; k - l)^,{x) - [(7."(s")')' - X'A V^^(^fc) - 9%{0) (2.5) 
r J 



i = 1, . . . ,N. Homogeneous system corresponding to (12.51) is 

y^ / aA;(a;)^j(a;)^.(x) - bk{x)ilj'j{x)ilJi{x) - Cfe(x)^j-(x)V^i(x) 



j=i 



+ — '?/'j(x)'?/'i(x) dx dj = 0, i = 1,2, . . . ,N 



(2.6) 



Let us multiply each equation in (12. 6p by d^ and add with respect to i: 

O'kix) I — ; ) — bk{x) — ; U]\[{x)dx — Ck{x)u1^{x) H — M^(x)) dx = (2.7) 



dx 



dx 



T 



As before, from (12. 7p it follows that un = 0, and therefore the homogeneous system (12.61) 
has only the trivial solution. This proves the uniqueness of the approximate solution un{x). 
Let us now prove uniform estimation of the sequence {un{x)}. Multiply (12. 4 p by di and add 
with respect to i = 1, . . . , A^: 



Sk 



0'k[X) [ —p-J - h{x)—^UN{x) - Ck{x)u%{x) + -ulf{x) + fk{x)UN{x) 



dx 



1 



Sk 

u{x- k - 1)um{x) dx - [(7s"(s")T - x]r] njv(sfc) + ^''ujv(O). (2. 



We estimate the four integrals on the left-hand side of (12. 8p as we did before to prove (12.30 
and derive 



dx 



1 Z"**^ 
dx ^ / u)^{x)dx <\g^\\uN{Q)\ + 
^^ Jo 



ls-{s-)'f\ + \xU]\uN{Sk)\ + 



Sk 



\fk{x)\ + -\u{x;k-l)\ 



T 



\u]^{x)\dx (2.9) 



for all r < ^. By Morrey's inequality we have 



max{|MAr(0)|; \uN{Sk)\} < ||M7v||c[0,.fc] < C\\uN\\wl[Q,Sk], 



(2.10) 



where the constant C is independent of N and r. By using Cauchy inequahties with appro- 
priately chosen e > 0, from (12.91) and (12.101) it easily follows that 

ii«ivirv^i[o,.,] < c(\\u{x- k - i)iu,[o,.,] + iiMu.(o,..) + i(7s"(s")')T + \xU' + i/r) (2.11) 

where C does not depend on A^, but depends on the time step r. From (12.111) it follows 
that {ujy} is weakly compact in W2^[0,Sfc]. Let v{x) be its weak limit point in iy2M0''^fc]- 
Passing to the limit in (12. 4p it follows that v{x) satisfies (I1.13P for ri{x) = ipi{x). Since 
{ipi} is a fundamental system in 14^2^[0,Sfc], it follows that v{x) satisfies (11.131) for every 
ri{x) G VF2^[0, Sk]. Hence v{x) = u{x; k) is a solution of (ll.lSp and in view of uniqueness the 
whole sequence un converges weakly in M^2^[0, Sk] to u{x; k). Lemma is proved. 
The following known criteria will be used in the proof of Theorem II. 2[ 

Lemma 2.2 IW^ Sequence of discrete optimal control problems In approximates the contin- 
uous optimal control problem I if and only if the following conditions are satisfied: 

(1) for arbitrary sufficiently small e > there exists number Ni = Ni{e) such that Qn{v) G 
V^ for all V G Vr^^ and N > Ni; and for any fixed e > and for all v G Vr-e the 
following inequality is satisfied: 

limsup (iNiQNiv)) - Jiv)) < 0. (2.12) 

7V-J.OO ^ ^ 

(2) for arbitrary sufficiently small e > there exists number N2 = A^2(e) such thatVN{[v]N) ^ 
Vfj+e for all [v]iy G V^ and N > N2; and for all [v]n G V^ , N > 1 the following in- 
equality is satisfied: 



limsup (j{Vn{[v]n))-In{[v]n)) < 0. (2.13) 



7V-s>oo 

(3) the following inequalities are satisfied: 

limsup j;(e) > J*, liminf j;(-e) < J*, (2.14) 

where J7'*(±e) = inf Jiu). 

Vr±, 

In the next lemma we prove that the mappings Q„ and Vn introduced in Section 11.31 satisfy 
the conditions of Lemma 12.21 

Lemma 2.3 For arbitrary sufficiently small e > there exists n^ such that 

Qniv) G V^, for all v G Vr_, and n > n,. (2.15) 

Vni[v]n) e Vr+„ for all [v]n G V^ and n > n,. (2.16) 

Proof. Let < e << R, v E V/j-e and Q{v) = [v]n = ([s]„, [g]n)- By applying Cauchy- 
Bunyakovski-Schwarz (CBS) inequality and Fubini's theorem we have 

n-l n-1 *'=+-'- X 

fc=i fe=i / ■^'^ 

tfc 

T t T T 

<lldtl W'iOl'd^ < I \s"{t)\'dt, J2^sl < I \s'{t)\'dt, (2.17) 



r 



t-T ''-^ 



*fc+l tfe 



n-1 „T n-l "J^ y 

Y.'^si- s\t)dt\ = \Y, / {s'ioyd^dt 

fc=o "^^ fc=o / •; 



< 



tfe+l i 



T 



'(■ -1- /* /* /* 

E / / [^'(^) + (^'(^))' ]^^^^ ^ ^ / [^'(^) + (^'W)' ]^^ ^ ^'^' (2.18) 

jfc=o r ; n 



ift ifc 



n „ n— 1 „2^ 

fc=l n fc=0 "^^ 



<i?V. 



(2.19) 



From ( l2Tri) -( l239|) it follows that 

max (\\ [s]n||^2, II [^]n||^) < max (||s||^2[o^r], ||^||^|[o,t]) + ^^^ < (^ - ^f + ^^^ < ^' 



2 

(2.20) 



if n > rZf 



RT 



1. Hence, (I2.15P is proved. 

Let us know choose [f]„, G V^. We simplify the notation and assume v = (s, g) = P„([f]n). 
Through direct calculations we derive 



IV1/|[0,T] 



n— 1 n— 1 n— 1 ^ ^ 



(2.21) 



fc=0 fc=l A:=l 

where C is independent of r. By using CBS inequality we have 

T T t 



TS- 



, < j \s'{t)\^dt, ^-sl, = L\J Js"{Od^dt 



< 







T t 



^J J\s"iO\'d^dt<^-J\s"it)\'dt. 



Since [v]n G V]^, from (E2ID,(E23) it follows that 
where Ci is independent of r. This implies that 



lim ||s| 



M^|[0,- 



n-l 



In a similar way we calculate 

k=0 k=l 

Hence, from (|2:2T]) . ( l2:22|) and fl2:25|) it follows that 

max (||s||^no^^], ||^||^^i[o,T]) < max (||[s]„||^2, ||M„||^) + Cr + - 

< R2 I ^ I "'"||q'||2 

^ K +L.T+ -\\S \\w^10,t] 



(2.22) 

(2.23) 
(2.24) 

(2.25) 



1,, n 



P Wwiio,- 



(2.26) 
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From fl2.24p it follows that given e > we can choose n^ such that for any n > Ue 

R' + Cr + ^\\sTw^^,^r] < {R + ef- (2.27) 

From (I2.26P and f l2.27p . ( I2.16P follows. Lemma is proved. 
Corollary 2.1 Let either [v\n € V^ or [t>]„ = Q„(f ) for t> G Vr. Then 

\sk - Sk~i\ < Ct, k = l,2,---,n (2.28) 

where C is independent of n. 
Indeed, if f G Vr, then s' G 14^2^0'^] ^^^ by Morrey inequality 

||s'||c[o,T] < C'lP'llvyjMo.T] - ^^R (2.29) 

and hence for the first component [s]„ of [v]n = Qn{v) we have (I2.28p . Also, if [v]n G V^, 
then the sequence f" = Vn{[v]n) belongs to Vr+i by Lemma [231 and the component s" of f" 
satisfies (I2.29p . Since, 

s"(0) = so, s^itk) = £^±f^, k = l,---,n 

from (|2:29D . (|2:28|) easily follows. 

3 Proofs of the Main Results 

3.1 First Energy Estimate and its Consequences 

Throughout this section we assume that 

4>eL2[0,l], j,xeW^^\D), aeL^{D), 

a satisfies (11.61) and b, c, / satisfy the conditions imposed in Section 11.41 The main goal of 
this section to prove the following energy estimation for the discrete state vector. 

Theorem 3.1 For all sufficiently small r discrete state vector [u{[v]n)]n satisfies the follow- 
ing stability estimations: 

'"' I du{x; k) 



max u'^(x;k)dx + Ty / , 
o<fc<n7o Z-^J^ \ dx 



dx < 



C[mU,so) + \\9\\U,T) + WfWliD) + ll7(^"(t),t)(^")'(t)llL(0,T) 

n—1 

+ llx(s"(t),t)||i,(o,T) + $^l+(^fc+i-5fc) r^\\x-k)dx), (3.1) 

fc=l "^^fe 

/"' 2/ , N , v^ /"' I duix] A;) 2 " /■' 

max / M (x: /e) da; + r > / ; dx + r y / Ujlx: kjdx < 

^^'^^-Jo t^Jo I dx f-io *' ' " 



fc=0''^ k=l 

C(\\<t>rwli0,s,) + ll^llL(o.T) + II/IIL(D) + ll7(^'^(t),t)(s")'(t)||i^(o, 

n-1 



T) 



+ llx(s"(i),t)llL(o,T) + I]l+(5fc+i-5fc) r^\\x-k)dx), (3.2) 

fc=i "^''fe 

where C is independent of r and 1+ 6e an indicator function of the positive semiaxis. 
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We split the proof into two Lemmas. 



Lemma 3.1 For all sufficiently small r, discrete state vector [u{[v]n)]n satisfies the following 
estimation: 



max / u^ix^k) dx + T} 



l<k<n 



fe=i 



Sk 



du{x; k) 



dx 



dx + T^y^ / Uj{x;k)dx< 



k=l 



CmUo,so) + ll^llw) + ll/lli.(D) + ll7(^"(t),t)(.7Wlli.(o. 



T) 



n-1 






(3.3) 



where C is independent of r. 

Proof. By choosing ri{x) = 2tu{x; k) in flL13p and by using the equahty 



2tu-[{x; k)u{x; k) = u (x; k) — u (x; fc — 1) + r %(a;; k) 



we have 



Sfc 



2 / „,2 
t 
JO JO JO 

k 



u{x;k)dx— / u {x; k — l)dx + T / Uj(x;k)dx + 2T ctfc(x) 



Sfe 



(iu(a;; /c) 



dx 



rfx 



2r / bkix 



du{x; k) 



dx 



u{x; k) + Ck{x)u (x; k) — fk{x)u{x; k) 



dx- 



2r [(7s 



n\l\k ^k 



X^n]u(sfc;A;)-2r/M(0;fc). 



(3.4) 



Using (II. 6p . Cauchy inequahties with appropriately chosen e > 0, and Morrey inequality 
( 1230|) from (jS^D we derive that 



Sk rsk 

u^{x]k)dx — I u^{x]k — \)dx + qqt 











ClT 



|(7..(5")ri' + lx."r + Kr+ / fk{x)dx+ u\x-k)dx 

Jo 



"^ du{x;k) 2 2 P 2/ ,^J / 

— dx + T / Uj[x; k)dx < 

ax Jq 

Sk PSk 



(3.5) 



where Ci is independent of r. Assuming thatr < Ci, from (13. 5p it follows that 



/"Sk 

(1 - Cir) / u^{x; k)dx < 
By induction we have 

••Sk 



Sk-1 



u^{x]k — l)dx + l^{sk — Sk-i) I u^{x;k — l)dx+ 







(3.6) 



psk fso "■ , 

/ u^{x-k)dx<{l-CiT)-^ (t)'^{x)dx + Y{l-CiT)-^^^-^\Ci 
Jo Jo j^i *- 



rx 



\{lsAsnyr + \xU' + \9'r+ f {x)dx +I4s,~s,^,) u'{x;j-l)dx\. (3.7) 



Sj-l 
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For arbitrary 1 < j < A; < n we have 

(1 - Cir)-'=+^-i < (1 - Cir)-'^ < (1 - Cir)-" = (l - ^ j ^ e^^^, (3.8) 

as r — > 0. Accordingly for sufficiently small r we have 

(1 - Cit)-''+^-^ < 2e^'^ ioTl<j<k<n, (3.9) 

By applying CBS inequality from (13. 7p - (13. 91) it follows that 

max ru\x-k)dx < C^lUWl,^,,,,) + 1|<7||L(o,t) + Il7(5"(t), t)(5")'(t)ilL(o,r) + 

Ms'\t),t)\\l(o,T) + \\f\\h(D) + J2^+isk+i-Sk) [ '^\\x-k)dx). (3.10) 

fc=l ^^k 

where C2 is independent of r. Having (13.101) . we perform summation of (13.51) with respect 
to k from 1 to ra and derive 

T" 2/ X , V^ P du(x;k)\ 2V^ P 2/ ,x, 

/ u [x; n) ax + aQT y^ — 1 ax + r /^ / Uj[x;k)dx< 

Jo j.^1 Jo ^ f^^i Jo 

ii0iiL(o,.o) +^i(ii^iiL(o,T) + wfwiiD) + Ms-it),t)isn'mi^o,T) 

" rsk ^ "'~^ rsk+1 

+ llx(s"(t),t)||i,(o,T)+^$^ / u\x;k)dx]+J2^+i-'^k+i-Sk) u\x-k)dx, (3.11) 

From (I3.10p and (13. lip . (13. 3p follows. Lemma is proved. 

In the next lemma we prove a nice property of the extension introduced in the Definition 
1.3, which allows to extend the estimation (13.31) to (13.11) and (13. 2p . 

Lemma 3.2 Given discrete control vector [f]„ G V^, a discrete state vector [u{[v]n)]n satis- 
fies the inequalty 



/ 2/ , X , V^ / I du(x; /c) 2 V^ r 2/ , X , 

max / u (x; k) dx + T y / dx + t y / u-r{x\k)dx< 

i<k<njr. '^-^ Jn \ dx ^^7n ~ 

/ r^ v^ r*' I rfwfx- A;) 2 J^ r'^k \ 

Ci max / u^{x]k)dx + Ty / — V^^ dx + r > / u|(a;; fc)cia; , (3.12) 

where C is independent of t. 

Proof. By induction it follows that the ffist two terms on the left hand side are estimated 
by the ffist two terms on the right hand side with the constant C = 2^, where A^ is defined 
in ( ITTiD . 

Define a family of functions {u{y; k),k = 0, ..., n} as 

m(2/; 0) = 0(i/so), uiy; k) = u{ysk; k), <y <l,k = l,...,n. 
As before, assume they are all continued by induction to semiaxis {y > 0} as 

u{y;k) = u{T-y-k), for 2""^ < y < 2". 
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We have 



n „ n „ 

^r / u^{x;k)dx<J2^ J 

A;=l n k=l n 



Uj{x; k)dx 



fc=i 



u{x/sk] k) - u{x/sk-i; k -1) 



T 



dx 



E 

k=l 



rsk 



u{y; k) - u{ysk/sk^i] k - I) 



dy<h + h 



(3.13) 



where 



/.=2E 



TSfe 



/c=l 



u{y]k) -u{y;k-l) 



r 



dy 



2^+^5^rs, ul{y;k)dy 



fc=i 



2JV+i^, 



fc=i 



u(x; k) — u{xsk-i/ Sk, k — 1) 



dx < 



Sk 



Sk 



k=l 



2^+2 J2 T / ulix; k)dx + 2^+2 ^ ^ 



k=l 



u{x; k — 1) — u{xsk-i/sk] k — 1) 



1 2 



dx (3.14) 



/. = 2E 



rsk 



k=l 



u{y; k-1) - u{ysk/sk-i; k - I) 



T 



dy. 



By using CBS inequahty, Fubini's theorem and Corollary 12.11 we have 



(3.15) 



Sk 



E 

k=l 



u{x; k -1) - u{xsk-i/sk] k-l) 



T 



dx 



Sk X 



E 

k=l 



du{C,; k — 1] 



■d^ 



Cfl 



n— 1 



d-^^H 



j. °fc-i 



fc=0 



(iM(x; fc) 



(ia; 



dx. 



(3.16) 



Ar2'' 
227V+1^2^2 - . 

fc=i •;; 



rfM(x; k-l) 2 23^+1^2 A^^ "4 A- ^^^. ^) 



(ix 



dx 



E 

fc=0 



dx 



dx 



23^+lC'2]V3/ ^ 

^ — A^E- 



fc=0 



^^(x; k) 



dx 



dx 



(3.17) 
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Hence, from (IXTSjl - dSTTj) it follows that 

I _, s 

Y^T fulix; k)dx < g(^r f\ ^''^^'' ^^ \dx + ^r j u^ix- k)dx^ (3.18) 



I 1 ^k ^k 



k=l g A:=0 Q fc=l 



where C is independent of r. From fl3.18p . p.l2p follows. Lemma is proved. 

It can be easily seen that Theorem 13.11 follows from Lemma 13.11 and Lemma 13.21 
Let [v\n G V'l^^n = 1,2, ... be a sequence of discrete controls. From Lemma [2.31 it follows 
that the sequence {"Pnliw]™)} is weakly precompact in iy2^[0,T] x iy2^[0,T]. Assume that 
the whole sequence converges to f = {s,g) weakly in 1^2^ [0,T] x iy2^[0,T]. This implies the 
strong convegence in 14^2^ [0,T] x L2[0,T]. Conversely, given control v = {s,g) G V^ we can 
choose sequence of discrete controls [f]„ = Q„(f). Appplying Lemma [53] twice one can easily 
establish that the sequence {PnlHn} converges to u = {s,g) weakly in iy|[0,T] x PF2^[0,T], 
and strongly in W^2^[0, T] x L2[0, T]. In the next theorem we prove the continuous dependence 
of the family of interpolarions {u^} on this convergence. 

Theorem 3.2 Let [v]n G V^, n = 1,2, ... be a sequence of discrete controls and the sequence 
{'Pn{[v]n} converges strongly in W}[0,T] x L2[0,T] to v = {s,g). Then the sequence {m'^} 
converges as r — t- weakly in W^ (fi) to weak solution m G V2 ' (fi) of the problem ^1.1\) - 
( |j.^| j, i.e. to the solution of the integral identity U.10\) . Moreover, u satisfies the energy 
estimate 

Ml^'\D) ^ C'(||0||i,(o,.o) + II^IIL(0,T) + II/IIL(D) + hWwl^^D) + \\^\\w^,'\d) (^'IS) 

Proof. In addition to quadratic interpolation of [s]„ from f ll.lip . consider two linear 
interpolations: 

~s^{t) = Sk-i + ^' ~ '""'' (t - tfe_i), tfc_i <t<tk,k = i;^; S"(t) = Sn, t> T; 
r 

5^(t) = 5'^(t + r), 0<t<r. 

It can be easily proved that both sequences S" and s" are equivalent to the sequence s" in 
14^2^ [0,T] and converge to s strongly in 14^2^ [0,T]. In particular, 

s^v\\si\\wi[o,T] < C* (3.20) 

n 

where C^, is independent of n. We estimate the last term on the right-hand side of (13. ip as 
follows: 

"~j- fSk + l "~^ ftk + l 

J2 u{sk+i - Sk) / u'ix- k)dx = Y^ i+(sfc+i - sk) j {r)'{t)u^{r{t)- k)dt = 

k=l "^^fe fc=l '^*fc 

P^k-\-~\ 

YUisk+i-Sk) / (n'W(«^(5"W,i-r)) dt = 
fc=i "^*fc 

YUisk+i-Sk) / (s7)'(t)(«-(s7(t),t)) dt. (3.21) 

1 — 1 'Jth~^ 



k=l ""^'=-1 
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By applying CBS inequality we have 



\Y,U{sk+i-Sk) f''^\\x;k)dx < ||(5D'IU.[o,T]||u^(5r(t),t)|lL[o,T]- (3-22) 

fc=i -^^h 

From the results on traces of the elements of space V2{D) ([23l HI [27] ) it follows that for 
arbitrary u G V2{D) the following inequality is valid 

\HrsUt),t)\\L4o,T]<C\\u\\v,iD), (3.23) 

with the constant C being independent of u as well as n. From (13.201) . (13. 221) and (13.231) it 
follows that 

\j2U{sk+i - Sk) u\x;k)dx <aC\\u^\\l^^j,y (3.24) 

fc=l -^^t: 

If the constant C* from (13.201) satisfies the condition 

C, < {CC)-^ (3.25) 

then from ([SH]) and iK2^ it follows that 



2 4- 

L2{D) + 



||7(3"(t),t)(.")'(t)||i^(o,^) + ||x(^"(t),t)||L(o,T)), (3.26) 



where C is another constant independent of n. By applying the results on the traces of 
elements of W2' (D) ([H [27]) on smooth curve x = s"(t), Morrey inequality for (s")' and 
([2A6|) we have 

||7(s"(t),t)(s")'(t)lU,(o,T) < \\{sn'\\cio,T]h{s''{t),t)\\L,io,T] < CMlw^.o^n) 

llx(^"(t),i)l|L.[o,T]<C3||x||^.0(^), (3.27) 

where C3 is independent of 7, x and n. Hence, from (13.261) and (I3.27P it follows the estimation 

lk1|J.x.0(^) < C(||0||i,(o,,„) + ||^||L(o,T) + II/IIL(D) + hWl^^o^o) + Mw^'\D))^ (3-28) 

with C being independent of n. 

If (I3.25P is not satisfied, then we can partition [0, T] into finitely many segments [tn _i, ^n ], 
j = l,q with t„o = 0, tng = T in such a way that by replacing [0,T] with any of the 
subsegments [t„ _-^,t„.] (13.201) will be satisfied with C* small enough to obey (13.251) . Hence, 
we divide D into finitely many subsets 

D' =Dr\{tn^_, <t<t„J 

such that every norm ||M'^||y c^,^) is uniformly bounded through the right-hand side of (I3.28p . 
Summation with j = 1, . . . ,q implies (I3.28p . 

From (I3.28P it follows that the sequence {u'^} is weakly precompact in W2' (-D). Let 
u G W2' (-D) be a weak limit point of u'^ in W2' (-D), and assume that whole sequence {m^} 
converges to u weakly in W2' (D). Let us prove that in fact u satisfies the integral identity 
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fll.lOP for arbitrary test function $ G W2 [fl) such that ^l^^j. = 0. Due to density of C^{U) 
in W2' (^) it is enough to assume $ G C^{Q). Without loss of generahty we can also assume 
that $ G C^(DT+r), ^ = 0, for T < t < T + r, where 

DT+r = {{x,t) : < X < /, < t < T + r} 

Otherwise, we can continue $ to D^+r with the described properties. 
Let 

$(x; A;) = $(x, fcr), $t(x; /c) = ^ ^ 

r 

As before, we construct piecewise constant interpolations $'^, $[. Obviously, the sequences 
{<l'^}, {^r} and {$[} converge as r — t- uniformly in D to $, |^ and ^ respectively. 

By choosing in fll.lSp ri{x) = r$(x;fc), after summation with respect to fc = l,n and 
transformation of the time difference term as follows 

n-l 



r ^ / Ut{x;k)^{x;k)dx = —Ty^ / u{x;k)^tix;k)dx— / 0(a;)$(a;; 1) rfx 
fc=l "^0 fc=l ■ " 

— y / m(x; fc)<l>(x; k) dx = — y 

k=l '' ^k k=l 

"~1 rtk+i pi rsyt) 

-J2 {~s'')'{t)u^{r{t),t-T)^^{s'^{t),t-T)dt=- / U^^dxdt 

k=i Jtk JO Jo 



tk~iJO 



u'^^l dx dt 



a;)$'^(x, r) dx 



rT ps{t) 



^0 



k r^k + 1 



fc=l Jtk-1 Js{t) 



M^$[ (ia; (it 



(3.29) 



we derive that 



^ Z"'"^*^ r 9m ■^ 9$^ (9m^ 1 f''° 

/ <^a- b^—<^^-cu^<^^ + f<^^-u^<^]}dxdt- (P(x)<^Ux,T)dx 

Jo [ dx dx dx J Jo 

\s'iy{tX{Cs^,){t),t)^mr,){t),t)dt+ [ gm^{o,t)dt 



+ 



7(s"(t), t)(s")'(t) - x(^"(t), i)) $^(s"(t), t) cit - i? = 



(3.30) 



where 



a— ^~^^ ~ CM *& + /$ > dxdt- y 






(9x 9x (9x 



fc r^k+i 



u'^^l dx dt 






a$" 



1 



7(s"(t),t)(s")'(t)-x(s"(t),t)) ^—dxdt+ / 0(x)$"(x,r)da; 



- 5a: 



so 
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Let 

n 

^ = U "l(^'^) • '^'t-i <t <tk, min(s(t),Sfc) < a; < max(s(t), s^)} 
fc=i 

I A| denotes the Lebesgue measure of A. Since 

" rtk rtk 2\/T 
|A| < ^ / / \s'{T)\dTdt< ||s'||l2(o,t)^^0 asr^O 

and all of the integrands are uniformly bounded in L^{D), it follows that the first term in the 
expression of R converges to zero as r — )■ 0. In a similar way one can see that the second and 
third terms also converge to zero as r — )■ 0. The last term in the expression of R converges 
to zero due to Corollary 12.11 and uniform convergence of {$^} in D. Hence, we have 

limi? = (3.31) 

Due to weak convergence of u'^ to u in W2' (D) and uniform convergence of the sequences 
1$'^}, {^r} and {$[} to $, || and ^ respectively, passing to limit as r — > 0, it follows 
that first, second and fourth integrals on the left-hand side of f l3.30p converge to similar 
integrals with u'^, $^, $[, $^(a;,r) and $^(0,t) replaced by m,$, ^, <l>(x,0) and $(0,t) 
respectively. Since s"" converges to s strongly in ^^2^0;^]; ^^^ traces 7(s"(t), (t)), x(s"(t),t) 
converge strongly in L2[0,T] to traces 7(s(t), (t)), x{s(t),t) respectively. Since $'^(s"(t),t) 
converge uniformly on [0,T] to $(s(t),t), passing to the limit as r — )■ 0, the last integral on 
the left-hand side of fl3.30p converge to similar integral with s^ and ^"^ replaced by s and $. 
It only remains to prove that 

/■T-T pT 

lim / {s1)\t)u-'{s1{t),t)^^{s1{t),t)dt= s\t)u{s{t),t)^{s{t),t) (3.32) 

Since {s"} converges to s strongly in PV2^[0,T], from 03.28P it follows that {n^(s"(t), t)} is 
uniformly bounded in L2[0,T] and 

||u^(s5'(t),t) -M^(s(t),t)||L2[o,T] ^0 asr^O (3.33) 

Since {u^} converges to u weakly in W2' (-D), it follows that 

M^(s(t), t) -^ u{s{t),t), weakly in LsiO, T] (3.34) 

Since {$^(s^(t), t)} converges to $(s(t), t) uniformly in [0, T], from flX^ .f lXMj) . flS:^ easily 
follows. 

Passing to the limit as r — > 0, from (l3.3Up it follows that u satisfies integral identity (ll.lOp . 
i.e. it is a weak solution of the problem (ll.ip - (ll.4p . Since this solution is unique (|23]) it 
follows that indeed the whole sequence {m^} converges to m G V2 ' (fi) weakly in W2' (fi). 
From the property of weak convergence and (I3.28p . (l3.19p follows. Lemma is proved. 

In particular. Theorem 13.21 implies the following well-known existence result ([23]): 

Corollary 3.1 For arbitrary v = {s,g) G Vr there exists a weak solution u G V^ ' (fi) of the 
problem U.1\) - [T^\ ) which satisfy the energy estimate Ii3.19\) 
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3.2 Second Energy Estimate and its Consequences 

The main goal of this section to prove the following energy estimation for the discrete state 
vector. 

Theorem 3.3 For all sufficiently small r discrete state vector [u{[v\n)]n satisfies the follow- 
ing stability estimation: 



max 

l<A:<n 



du{x; k) 



dx 



n pi 




fc=l 



dx + T yj / u^{x; k) dx < 



C 



^.MM + ll^li;|,,,, + ll7(^"W,^)(^") ' 



+ llx(s"(t),t)ir 1 + 



n-l 



+ 



+Ei- 



L2iD) -r /_^ ^+[Sk+i - Sk) I u^ix; k)dx 
k=l -^^k 



Sfc + 1 



(3.35) 



We split the proof into two lemmas. 

Lemma 3.3 Let given discrete control vector [v\n, along with discrete state vector [u([f ]„)]„,, 
the vector function 

[u{[v]n)]n = iu{x; 0), u{x; 1), ..., u{x] n)) 

is defined as 

~, . u) _ i "(^; ^) < X < Sfc, 

\ u{sk]k) Sk < X < l,k = 0,n. 

Then for all sufficiently small r, [u{[v]n)]n satisfies the following estimation: 



max 

l<fc<n 



c 



du{x; k) 



dx 



k=l 



dx + t\^ / u\{x]k)dx + t'^2. 



fc=i 



Sk 



du{x; k) 
dx 



dx < 



w^oi\ + \\9f \ +ll7(s"(t),t)(s7Wll' i +llx(s"(t),t)f 1 + 



+ 



L2{D) 



^l+(sfc+i -Sfc) / u^{x]k)dx 

1 1 J Sl. 



r^k+i 

k=l "'^fc 

Proof: By choosing ?7(x) = 2r-uj(x; k) in fll.l3p and by using the following identity 



(3.36) 



, ,du(x;k) /du(x;k)\ , , f du{x]k)\ , , / du(x; k — 1)^ 

2rafc(x) 1 I = ak{x) { ) - ak-i{x} 



dx 



dx 



-TOkiix) 



du{x] k — V 
dx 



dx 



+ r^afe(a;) 



dx 



du{x; k) 
dx 



(3.37) 
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we have 



du(x;k)\ 
o-k[x] I -, I dx 







dx 

+T^ / afc(x) 
Jo 

+2t / hix) 
Jo 

Sk 



ak-i{x) 



du{x; k — 1) 
dx 



du{x; k) 
dx 
du{x; k) _ 



Sk 

<T I akt{x) 



dx + 2t / {utix; k)) dx+ 
Jo 

2 



du{x; k — 1] 
dx 



dx 



dx 



Ut{x; k) dx + 2t / Ck{x)u{x; k)ut{x; k) dx 



-2r / Mx)ut{x; k) dx - 2r [(%"(s")')' - x'A nisk] k) - 2r(7fc%(0; k) (3.38) 



By adding inequalities fl3.38p with respect to k from 1 to arbitrary m < n we derive 



(^mX,-^ ) 



< T 



du{x; m) 
dx 

"^ rsk 
X^ / afct( 



k=l 



r^k 

dx + 2r y^ / Mf (x; A;) dx + r^ y 



k=l 



ak[x) 



du{x; k) 
dx 



tJ 



dx 



k=l 







du{x; k — 1) 



dx 



dx + 2t\^ / bk{x) — j-^ — Ut{x;k)dx 



k=l 



dx 



+2tJ2 



fc=l 



Sk 



Ck(X uix 



'"■ PSk 

; k)ut{x; k) dx — 2r N^ / fk{x)ut{x] k) dx 



fc=i 



so 



+ j «o(a;)l^l rfx-2r^[(7,n(s"y)'=-xy iif(sfe;fc)-2rJ](7fcMf(0;A;) (3.39) 

By using (ll.6p .( irT6|) and by applying Cauchy inequalities with appropriately chosen e > 0, 
from (13.391) it follows that 



oo 



du{x; k) 



dx 



dx + T y^ / Uf-{x; k) dx + a^r"^ 2. 



k=l 



k=l 



Sk 



du{x; k) 
dx 



dx < 



m r 

<ctE 

fc=l 



u^{x] k) dx + 



Sk 



du{x; k) 



dx 



Sk 



dx + fj^ (x) dx 



C f- dx-2TY,hs-{s^)'f-x':A^t{sk\k)-2rY,9kUt{'^\k)dx (3.40) 

"^o ^ fe=i fc=i 

where C is independent of n. Note that we replaced u with u in first two integrals on the 

right-hand side of ^M^. Since 7, x e Wl'^{D) we have 7(s"(t), t), x(s"(t), t) G 1^2^ [0,T] 
([211111 [23]) and 

where C is independent of n. According to Lemma I2T31 P.. ( [f ] ^ ) G Vr+i. By applying Morrey 

1 

inequality to (s")' we easily deduce that 7(s"(t), t)(s"')'(t) G 1^2* [0;^] ^^^ moreover, 



Us-{t),t){s^ 



w^ [o,r] 



<Ci||7(s"(t),t)|| 



W^ [0,T] 
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where C is independent of n. 

Let w{x,t) be a function in W2' (D) such that 

w{x, 0) = 0(x) for X G [0, sq], a(0, t)«;:,.(0, t) = c/(t), for a.e. t G [0, T] (3.43) 

a(s"(t),t)w,.(s"(t),t) = 7(s"(t),t)(s")'(t) - x{s''{t),t) for a.e. t G [0,T] (3.44) 



and 



|W||^2,1(^) < C 



+ ||7(."(t),t)(.")'(t)-x(^"(t),t)||^^.^^^^^ 



so] 



(3.45) 



The existence of w follows from the result on traces of Sobolev functions [U HTj. For example, 
w can be constructed as a solution from W2 ' {^^) of the heat equation in 

fi" = {0<x<s'^(t),0<t<T} 

under initial-boundary conditions (I3.43p .f l3^4^ with subsequent continuation to W2' (-D) with 
norm preservation [32l |33] . 

Hence, by replacing in the original problem (ll.ip - (ll.4p u with u—w we can derive modified 
(I3.40p without the last three terms on the right-hand side and with /, replaced by 



F = f + wt- (aw^)^ - bw^ -cw e L2{D). 



(3.46) 



By using the stability estimation (13. 3p . from modified (13. 40 p . ( 13^451) and (I3.46p . (I3.36P follows. 
Lemma is proved. 

In the next lemma we prove (I3.35P with I being replaced with Sk on the left-hand side. 

Lemma 3.4 For all sufficiently small t, discrete state vector [u{[v]n)]n satisfies the following 
estimation: 



max 

l<fc<n 



Sk 



du{x; k) 



dx 



rsk 
dx + ry^ / uf{x;k)dx< 



C 



fc=i 

„n\i / ±\ i|2 



wM^i] + \\9\? \ +Il7(s"(t),t)(s")'(t)|r 1 +||x(s"(t),t)f 1 + 



+ 



2 
L2{D) 



'"■ 1 rsk+i 

+ ^l + iSk+l-Sk) j 

fc=i "^''fc 



u^{x] k)dx 



(3.47) 



Proof: Obviously, we can equivalently replace u with u in the first term on the left-hand 
side of (I3.36p . We can do so also in the second term provided Sk-i > Sk for all k = l,m. 
Hence, we only need to estimate 



uj{x;k)dx, Sk-i < Sk- 
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By using (12.281) we have 



Sk 



«*-! Sk 



Uj{x;k)dx = / Ui[x]k)dx + / u^{x;k)dx 

Sk-i 



Sk 



u{x; k) — u{x; k — 1) 



Sk 



dx <2 



u{x; k) — u{2sk-i — x; k) 



dx 



Sk-l 



Sk 


u{2sk-i -x;k)- u{2sk-i - 

T 


X 


;fc-i) 


2 


Sk 

dx<2 


X 

1 /■ du{y; k) 
r J dy 


Sk-l Sk-l 2sk-i—x 


Sk-l Sk X 

+2 / M^(a;; k)dx <— / 


du{y; k) 
dy 


2 

dy2{x — Sk-i)dx 


Sk-i-{sk-Sk-i) Sk-l 2sk-i-x 


Sk-l Sfe 

+2 / ul{x; k)dx < 2 / 


du{x] k) 
dx 


Sk-l 

dx + 2 I u\{x] k)dx. 


Sk-i-Cr Sk-i-Cr Sk-i-Cr 


i, for sufficiently small r we have 




Sk Sk 

ul{x;k)dx<2 

Sk-i-C- 


r 


du{x\ k) 
dx 


2 
( 


Sk-l 

ix + / u^{x]k)dx 





dx 



(3.48) 



Sk 



+2 / ul{x] k)dx < 2 

Sk-i-Cr 



du{x; k) 



dx 



dx + ?> I u^{x; k)dx 



(3.49) 



From (I3.36P and (I3.49p . (I3.47P follows. Lemma is proved. 

It can be easily seen that Theorem 13.31 follows from Lemma 13.41 and extension Lemma 13.21 
Second energy estimate (I3.35P allows to strengthen the result of Theorem 13.21 

Theorem 3.4 Let [v]n € V^,n = 1, 2, ... be a sequence of discrete controls and the sequence 
{Vn{[v]n} converges strongly in W2[0,T] x L2[0,T] to v = {s,g). Then the sequence {u'^} 
converges as r — )■ weakly in W2 ' (fi) to weak solution u G W2 ' (fi) of the problem M-l\) - 
( [i.^| ), i.e. to the solution of the integral identity U.9\) . Moreover, u satisfies the energy 
estimate 



l7/l|2 < r 



2 I |U||2 I II /l|2 



W^[0,T] 



L2{D) + \h\\w^'\D) + ll^llwi'i 



1^2 •(^) 



(3.50) 



Proof: The last term on the right-hand side of the second energy estimate (I3.35P is esti- 
mated in Theorem O along fl3:20l) - fl3:24l) . By using Theorems O and 1121 from (JHJS]) , (JSJl 
it follows that the sequence {u^} satisfies the estimate 



\^^\\w^-\d) - ^ 



WiiO,so) + ll^llp^ip^T'] ^ "•^"^2(^) + ll^llw"m^ + 11X11 wM 



|2 



2 
WnD) 



(3.51) 



22 



Hence, {u'^} is weakly precompact in W2' [D). It follows that it is strongly precompact in 
L2{D). Let M be a weak limit point of {u'^} in W2' (D), and therefore a strong limit point 
in L2ID). From (K^ it follows that 

1 . % /■ 

\\u'^ — m'^II^^j-j^) = -r^ y^ / Mf (x; /i;)(ia; — )■ 0, as r — )■ 0. 

'==10 

Therefore, m is a strong limit point of the sequence {m^} in L2{D). By Theorem 13.21 whole 
sequence {u^} converges weakly in W2 {^) to the unique weak solution from V2 ' {^) of the 
problem fll.ip - fll.4l) . Hence, m is a weak solution of the problem fll.ll) - fll.4p and we conclude 
that whole sequence {u^} converges weakly in W2 (D) to m G W2' (D) which is a weak 
solution of the problem f ll.ip - fll.4l) from W2' (fi). From the property of weak convergence it 
follows that u satisfies (13.501) . Theorem is proved. 

In particular. Theorem 13.41 implies the following existence result: 

Corollary 3.2 For arbitrary f = (s, (?) G Vr there exists a weak solution u G W2' (fi) of the 
problem / li.i|) - PT^ which satisfy the energy estimate \3.50^) 

Remark: In fact, we proved slightly higher regularity of m, and both in Theorem 13.41 and 
Corollary 13.21 Wn ' (-D)-norm on the left-hand side of fl3.50p can be replaced with the norm 

v<t<l ^ 

3.3 Proof of Theorem [TTTl 

Let \yr^ G Vr be a minimizing sequence 

lim J{vn) = J* 

n—^oo 

Sequence w„ = (s„,(7„) is weakly precompact in iy2^[0,T] x 1^2^ [0,T]. Assume that the 
whole sequence f„ = (s„,f7„) converge to some limit function v = {s,g) G Vr weakly in 
W^[0,T] X W^[0,T], and hence, strongly in Wi[0,T] x L2[0,T]. Let m„ = u{x,t;Vn),u = 
u{x,t;v) G W2' (-D) are weak solutions of (ll.ip - flL4p in W2' (fin) and W2' (fi) respectively. 
By Corollary 13.21 both satisfy energy estimation fl3.50p with Qn and g on the right-hand 
side respectively. Since f„ G Vr, \\un\\-^Yi.irj^\ is uniformly bounded. Hence, the sequence 
Am = Un — u it satisfies 

||Am||^i,1(^) < C (3.52) 



"2 

-1,1/ 



uniformly with respect to n. Accordingly, {Am} is weakly precompact in W2' (D). Without 



1,1, 



loss of generality assume that the whole sequence Un—u converges weakly in W2 ' (D) to some 



function v G W2' (D). Let us subtract integral identities (11.90 for Un and u, by assuming 



rl,l 



that the fixed test function $ belongs to W2' (D). Indeed, otherwise $ can be continued to 
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D as an element of W^ [D). 

/ <^ aAn^$^. - 6Am^.<I> - cAu^ + Aut^^ >dxdt+ / (5f„(t) - g{t)) $(0, t) dt 

+ / [7(s„(t), t)<(t) - 7(s(t), t)s'(t) - x{sn{t),t) + x(s(t), t)] $(s(t), t) rft 
Jo 

+ / / {a{un)x'^x - b{un)x'^ - cn„<l> + (M„,)i$ + /$} da; dt 

Jo Js{t) 
+ [ [^{Sn{t),t)s'„it)-x{Sn{t),t)]mSn{t),t)-<!>{s{t),t)]dt = (3.53) 

By using energy estimate fl3.50p . and continuity of traces 'y{s(t),t),x{s(t),t) of elements 
7, X ^ W2' (-D), strongly in L2[0,T] with respect to s G iy2^[0,T], passing to the limit as 
n —7- +00, from fl3.53p it follows that the weak limit function v satisfies the integral identity 

/ / {avx^x - bvx^ + cv^ + vt^} dxdt = (3.54) 

Jo Jo 

for arbitrary $ G W2' (D). Since, any element $ G Wo' (fi) can be continued to D as 
element of W2' {D), (13.541) is valid for arbitrary $ G W2' {^)- Hence, f is a weak solution 
from 1^2 '^(fi) of the problem ([III])-([LlD with f = g = -f = x = 0. From fl330D and 
uniqueness it follows that v = 0. Thus u„ converges to u weakly in W2' (-D). From Sobolev 
trace theorem ([4ii27]|) it follows that 

||M„(0,t) -u{0,t)\\L'2[o^T] -^0, ||M„(s(t),t) -u{s{t),t))\\L'2[o,T] ^0 aS Tl -^ OO, 
\\UniSnit),t) -M(s(t),t))||i2[o,T] < ||M„(s„(t), t) " M„(s(t), t) ||l2[o,t] 

+ \\un{s{t),t) - u(s(t),t))||i2[o,r] -^0 as n -^ oo. 

Hence, we have 

J{v) = lim J{Vn) = J, 

n—>oo 

and V is a solution of the Problem I. Theorem is proved. 

Remark: By applying first and second energy estimates we proved that functional J'{v) 
is weakly continuous in 14^2^ [0,T] x VF2M0'^]- Since Vr is weakly compact existence of the 
optimal control follows from Weierstrass theorem in weak topology. 

3.4 Proof of Theorem [7721 

We split the remainder of the proof into three lemmas. 
Lemma 3.5 Let J7i,(±e) = inf J{v), e > 0. Then 

lim j;(e) = j; = lim j;(-e) (3.55) 

Proof: Note that for < ei < £2 we have 

J*{t2) < JM) <J*< M-ei) < j;(-e2) 
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Therefore lim jr*(e) < J^ and limj7*(— e) > J^ exist. Let us choose f^ G Vp+e such that 

e->-0 e-S>0 



lim(:r(tO- j;(e)) =0 (3.56) 

Since f^ = (s^, g^ is weakly precompact in Vr2^[0, T] x VF2^[0, T], there exists some subsequence 
e' such that 

s^t — !■ s^, weakly in W2^[0,T], g^i — )• g^ weakly in PFg^pjT] as e' — )• 

Moreover, we have v^ = {s^,g^) G Vr. Since J7'(f) is weakly continuous it follows that 

\imJM = J{v,). (3.57) 

e'->-0 

From f l336|) . ( 13371) it follows that 

lim j;(e') = J* 

which implies the first relation in 03.55p . 



To prove the second relation in (I3.55p . take eo > and v = {s,g) G Vn-eo- Let {ak} be a 

-, lim Ofc = and set 

k—^+oo 

Vk = (sk, gk) = akV + (1 - ak)v^ 



real sequence with < a^ < 1, lim a^ = and set 

fc— >+oo 



where J{y^) = J^. We have Vk G V^-akf-Q and tifc converges to f* strongly in 14^2^ [0,T] x 
14^2^ [0,T]. Since J{v) is continuous, Vk is a minimizing sequence: 

lim J{vk) = J* 

fc— >oo 

For fixed k choose arbitrary e such that < e < egafc. We obviously have 

J*{-e) < J{vk), < e < eoOk 
Passing to the limit as e — )■ we have 

limj;(-e) <J{vk) 

e->-0 

Now we pass to the limit as A; — > +oo and get 

limj;(-e) < j; 

e->0 

Since the opposite inequality is obvious, (I3.55P follows. Lemma is proved. 
Lemma 3.6 For arbitrary v = {s,g) G Vr, 

lim In{Qn{v)) = J{v) (3.58) 

n— >oo 

Proof: Let v G Vr, u = u{x,t;v), Qn{v) = [v]n and [M([f]n)]n be a corresponding discrete 
state vector. In Theorem 13.41 it is proved that the sequence {u^} converges to u weakly 
in W2' (fi). This implies that the sequences of traces {u'^{0,t)} and {u'^{s(t),t)} converge 
strongly in L2[0,T] to corresponding traces u{0,t) and u{s(t),t). Let us prove that that 
the sequences of traces {u'^{0,t)} and {u'^{s(t),t)} converge strongly in L^[0,T] to traces 

25 



u{0,t) and u{s{t),t) respectively. By Sobolev embedding theorem ([ll[2^) it is enough to 
prove that the sequences {u'^} and {it'^} are equivalent in strong topology of W2' (fi). In 
Theorem 13.41 it is proved that they are equivalent in strong topology of L2{D). It remains 
only to demonstrate that the sequences of derivatives {ul} and {itl} are equivalent in strong 
topology of L2(fi). We have 



^ mm(sfc_i;Sfc) 



I T "T||2 



4^^' f (^)'<'-+ll<-"^llMr.). (3-59) 

U 1 ^ 



k=l 



where Sk = s"'{tk), s" is the first component of Vn{[v]n) and 

r„ = Uk=i{tk-i <t<tk, min(sfc_i; Sk) < x < s{t)} 

Since s" converges to s uniformly on [0,T], it follows that the Lebesgue measure of r„ 
converges to zero as n — )■ +00. By Theorems 13.21 and 13.41 the integrand is uniformly 
bounded in L2{D). Therefore, the second term on the right-hand side of f l3.59p converges 
to zero as n — ?■ +00. First term on the right-hand side of fl3.59p converges to zero due to 
stability estimation fl3.36p and the claim is proved. 

Let jy'^{t) = ly'^, fJ-^it) = /i'^ , if 4-1 < t <tk, k = 1,. . . ,n. We have 

||z/'' - z/||j^2[o_r] -)■ 0, ll/i^ - /i||i2[o,r] -)■ as r -)■ (3.60) 

We estimate the first term in X„(Q„(t')) as follows 

/3ot^|m(0;A;)-z/Y = /3oX1 / HO;k) - u''^ dt = ^o |m"(0, t) - z/"(t)p rft (3.61) 

j = l k=l -^^k-l Jo 

From IK60\i it follows that 

n 

lim /3or V |m(0; k) - z/^ = /3o||w(0,t) - z/(t)l|i2[o,T] (3.62) 

fc=l 

We estimate the second term in X„(Q„(t>)) as follows 

/9ir^|M(sfc;A;)-/iY = 2/3i5^ / / ^- {u^{s{t),t) - fj^it)) dxdt 
k=i k=iJ^>=-iJs{t) ox 

+A Y / Hs{t)- k) - fi'^l^ dt + PiY / -T^dx] dt = h + h + h (3.63) 

t^lJik-i k^l-Jtk-1 \Js{t) ox J 

We have 

limJ2 = /3i/ \u^{s{t),t)-fi^{t)\'^dt = (3il \u{s{t),t) - n{t)\'^ dt (3.64) 



Since ||(w^)i-||l2(d) and \\u'^{s(t),t) — fi'^\\L2[o,T] are uniformly bounded, and {s"} converges 
to s uniformly on [0,T], by applying CBS inequality it easily follows that 

lim h = 0, lim J3 = (3.65) 

n—^oo n—^OD 
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From fl3:63|) -( l335|) it follows that 

n „T 

\iml3iT^\u{sk]k)- ii^\^ = (3i j \u{s{t),t)- n{t)f dt (3.66) 

"""^ fc=i -^0 

Therefore, from (13.611) and (I3.66p . (I3.58p follows. Lemma is proved. 
Lemma 3.7 For arbitrary [v]n G V^ 

hm (j{Vn{[v]n)) - X„(H„)) = (3.67) 

Proof: Let [f]„ G V^ and f" = (s",*?") = P„([f]„). From Lemma [231 it follows that the 
sequence {'Pn(b]n} is weakly precompact in iy2^[0,T] x 1^2^ [0,T]. Assume that the whole 
sequence converges to w = {s,g) weakly in iy2^[0,T] x H^2^[0,T]. This implies the strong 
convegence in 14^2^ [0,T] x L2[0,T]. From the well-known property of weak convergence it 
follows that V & Vr. In particular s" converges to s uniformly on [0, T] and we have 

lim max |s"(ti) - s{ti)\ = (3.68) 

n—^oo 0<i<n 

Let Qn{v) = [v]n We have 

In{[v]n)~J{vn=MM -Xn([^]n) + In{l^]n) - J{v) + J{v) - Jiv'^) (3.69) 

In Section 1X51 we proved the weak continuity of the functional J^{v), i.e. 

lim {J{v) - Jiv"")) = 0. 

n— >oo 

From Lemma 13.61 it follows that 

lim {I^{[d]n) - J{v)) =0. 

n— >oo 

Hence, we only need to prove that 

lim {Xn{[v]n)-ln{[v]n))=0 (3.70) 

n—>-oo 

Let 

[«(Hn)]n = {un{x; 0), M„(x; 1), ...,Un{x; u)) , [u{[v]n)]n = {u{x; 0),u{x; 1), ...,u{x; n)) 

are corresponding discrete state vectors according to Definition 11.31 Let Sk = s"(tfc), Sk = 
s(tfc) and 

Au{x; k) = Un{x; k) — u{x; k) 

We have 

n n 

X„(H„) -X„([{)]„) = /3o$^r (n„(0; fc) - /o')' + A J]r (n„(sfc; fc) - /f)' 

fc=i fc=i 

n n n 

-(],J2r{mk)-f^Y-PiJ2T{u{h;k)-f^Y = PoJ2r{Au{0;k)f + 

k=l k=l k=l 

n n 

+2/3o Y, ^^«(0; k) (^(0; k) - /o ) + /^i 5Z ^ {^u{su\ k) + u{sk; k) - u{sk\ k)f + 

k=l k=l 

n 

+2/3i Y^ T (Au(sfc; k) + u(sfc; A:) - u{~Sk] k)) {u{~Sk] k) - /f ) (3.71) 



fc=i 
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From the estimations of Sections 13.11 and 13.21 it follows that the sequences {un{x;k)}, 
{u{x; k)} are uniformly bounded in ^2^(0,/]. From fl3.68p it follows that 



^ T {u{sk] k) - u{sk; k)f = ^ 



r 



fe=i 



fc=i 



Sfc 



Sfc 



du{x] k) 
dx 



dx 



< max Isfc — Sfcl > r 



l<fe<n 



fc=l 



Sk 
Sk 



du{x; k) 



dx 



dx 



— )■ 0, as n — )• +00. 



(3.72) 



From f l3.7ip and f l3.72p it follows that in order to prove fl3.70p it is enough to prove that 



R = J2^\{Au{0;k)f+{Au{sk;k)y 



k=l 



— 7- as n — )■ +CXD 



(3.73) 



By the Morrey inequality we have 



i?<C^r 



fc=i 



Am(x; A;) dx 



Sfc 



dAu{x; k) 



dx 



dx 



(3.74) 



where C is independent of n. Let us subtract integral identities 01.13P for Un{x] k) and 
■u(x; fc), by assuming that the fixed test function t] belongs to W2^[0,/]. Indeed, otherwise i] 
can be continued to [0, 1] as a element of W^2^[0, /]: 



"^^ f dAu dri , , . dAu ,, ,,. . \, k t / \ /~\n 

I ^>'^^>~J~1 Ok[x)——ri[x) - Ck[x)Auri + AutT] \ dx - Xs msk) - ri[Sk)\ + 

+ / \(^k{x)--- bk{x)-—r] - Ck{x)uri + fk{x)ri + UtT] ] dx+ 

Js^ \ dx dx dx J 

+ [{is^is'^yf - {Tssf] Visk) + (Ts-sT ivisk) - Vih)] - [x'n - Xl] Visk) = (3.75) 



Our goal now is to derive from (I3.75P that the right-hand side of (I3.74p converges to zero as 
n — )■ +00. The proof goes along the same lines as the derivation of the first energy estimate in 
Lemma [3m By choosing ri{x) = 2tAu{x; k) in (13.751) . and by using (II. 6p . Cauchy inequalities 
with appropriately chosen e > 0, and Morrey inequality (I2.10p we derive similar to (13. 5p : 



QqT 



Sk 



dAu{x\ k) 



dx 



Sk 



Sk 



"Sk 

-2 / A„,2/ 



Sk 



dx + Au{x;k)dx— / Au {s;k — 1) dx + t / Au^{x;k)dx 

Jo Jo Jo 

Sk 



<CiT / Au\x]k)dx-2T 

Jo J Sk 



+fk{x)Au{x; k) + UtAu{x; k) 



dudAu(x;k) , , .du . . ,, / x~a / t\ 

(^kix)- bk{x)—Au{x; k) - Ck[x)uAu[x; k) 

dx-2 r {-f{s'\t),t){s'')\t)--f{5{t),t)S\t))dt Au{sk]k) 
Jtk-i 

-2 f" r\(~s(t),trs'{t)^^^^p^dxdt + 2 r {x{s''{t),t) - x{m,t)) dtAu{sk;k) 

Jtk-l Jsk "^ Jtk-l 



dx 

['" r" , , , ^dAu(x;k) , , 
+2/ / x(s(t),t) y^^dxdt 

Jtk-l Jsk 



dx 



(3.76) 
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By applying the technique along with f l3.6p - fl3.10p from fl3.76p it follows that for all sufficiently 
small T 

A-u^(x; A;) (ix < C( >^ l+(sj+i — Sj) / A'u^(x; j)(ix + >^ |£j| j (3.77) 

.7 = 1 ^^^ .7 = 1 



where C is independent of r and 

du dAu{x;j) 



Cj = T 



"^•^^^rfx dx 



dvL 
- bjix)—Au{x;j) - Cj{x)uAu{x;j) + fj{x)Au{x;j) 



+UtAu{x;j))dx+ / [7(s"(t),t)(0'W -7(5(t),t)5'(t)]cit Au{s,;j) 



*7-l 



Ij _ 1 ^ S j 



7(5(t),t)5'(t)^^^^^^dxdt+ /"' (x(s"(t),t)-x(5(t),t))dt Am(s,;j) 



dx 



•j-i 



+ 



j-i "^ *j 



X(s(t) , t) — ^|l^ c?a; c?t 



(3.78) 



Having fl3.77p we perform summation in fl3.76l) with respect to k from 1 to n and derive 



l<fc<n 



max / Au'^(x;k) dx +} r 

\<k<n L \ ' / /_^ 



fc=i 



dAu{x; k) 



dx 



n-l 



dx + y^ r^ / Auf-{x; k) dx 

k=l "^0 

<Ci(5^1+(s,+i-s,) / Au\x;j)dx + J2\^^\) (3-79) 

i=i "^^^ i=i 

where Ci is independent of r. Our next goal is to absorb the first term on the right-hand side 
of fl3.79p into the left-hand side. We apply the same method used in the proof of Theorem l3.2l 
(see fl3.20p - fl3.26p ). The only difference is that in the estimations fl3.23p and fl3.24p we replace 
D with 

Ql = {0 <t <r,0<x<s^(t)} 

Let us also introduce the region 

n 
fin = [J {4-1 <t<tk,0<X<Sk} 

Note that 



fc=i 



I Am" II 2 1,0 



Vj ' (^n) l<k<n 



max / Au^{x\k)dx+y ■ 

\<k<n L \ ' / /_^ 



fc=i 



Sfc 



dAu{x; k) 



dx 



dx 



Hence, by applying the method used in Theorem 13.21 we derive from fl3.79p the following 
estimate: 



j=i 



+ EIA 



(3.80) 
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where C2 is independent of r. Since the sequences {s"} and {s"} are equivalent in strong 
topology of VF2^[0,T], the first term on the right-hand side of fl3.80p converges to zero as 
n — )■ +00. It only remains to prove that 



lim y^\CA =0. 



n— >+oo 



(3.81) 



We have 



E 



i=i 



J du{x; j) dAu{x; j) 

T I a Ax) ; ; dx 

dx dx 



<C 



diT' 



dx 



L2iA) 



E 

dAu 



dx 



t,_i Js^ 



L2(A) 



du'^ dAu^ 
a{x, t)— — dx 



(3.82) 



where 



A = M{(a;, t) : tj-i <t< tj, min(sj, Sj) < x < max(sj, Sj)} 
3=1 

From (13.681) it follows that the Lebesgue measure of A converges to zero as n — )■ 00. Since 
by the first energy estimate W2 ' (-D) norm of u'^ and Am"^ are uniformly bounded, the right- 
hand side of (I3.82P converges to zero as n — )■ cxd. For the same reason, the next three terms 
in the expression of X]?=i I A' I ^^^^ converge to zero as n — >• 00. We have 



E 



r / Ut{x;j)Au{x;j)dx 



- (5Z^ / '^Kx;J)dx] II Am 



Il2(A) 



(3.83) 



By the second energy estimate W2 ' (-D) norm of u'^ is uniformly bounded. Accordingly, the 
right-hand side of (I3.83P converges to zero as n — )■ +00. We have 



n rt 



E, 



n „/ . 



7(s"(t),t)(s")'(t)-7(5(t),t)5'(t)jAM(s„j) 
7(s"(t), t)(s")'(t) - 7(S(t), t)5'(t)) f An(S(t); j) + 



dAu^{x,t) 



m 



dx 



dx 



dt 



<||7(s"(t),t)-7(5(t),t)|U,[o,T]||(s")1lc;[o,T]||An-(S(t),t)|U,[o,T] 
+ ||(.")'(t)-5'(t)|U,[o,T]||7(5(t),t)|U4o,T]||A«-(5(t),t)|U,[o,T] 



^ / rt 



E / l7(^"(t),t)(.")'W -7(5(t),t)5'(t)rrft 



■j-i 



s(t) 



(9An^ 
dx 



dx 



2\ 2 



< C||7(s"(t),t) -7(5(t),t)|U,[o,r]||(s7lk,Mo,T]l|A^^1 



W^2' («) 



+C||(.")'(t) - 5'(t)|U,[o,T]||7lU,^.v.JIA«^ 



+ ||7(s"(t),t)(.")'(t)-7(5(t),t)5'(t)|U2[o,T] 



i2[0,T]||7llT^2i.i(D)l 

9Am^ 



kaCI?) 



9x 



L2{A) 



max |5(t)-s"(t)| (3.84) 

0<t<T ' 
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and all three terms on the right-hand side converge to zero as n — )■ +00. Similarly one can 
prove that all the last three terms in the expression of X]?=i I A' I converges to zero as n — )■ cxd. 
Hence, fl3.8ip is proved. From f l3.80p and fl3.8ip . f l3.73p follows. Lemma is proved. 
Having Lemmas 13. 5[ 13.61 and 13. 7[ Theorem 11.21 follows from Lemma 12.21 
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